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ABSTRACT 

Aims. We use the Barnard 68 dark globule as a test case for a spherically symmetric PDR model exposed to low-UV radiation fields. With a 
roughly spherical morphology and an accurately determined density profile, Barnard 68 is ideal for this purpose. The processes governing the 
energy balance in the cloud surface are studied in detail. 

Methods. We compare the spherically symmetric PDR model by Storzer, Stutzki & Sternberg (1996) to observations of the three lowest 
rotational transitions of ''CO, '''CO / = 2 — > 1 and / = 3 — > 2 as well as the [Ci] 'Pi-^'Po fine structure transition. We study the role of 
Polycyclic Aromatic Hydrocarbons (PAHs) in the chemical network of the PDR model and consider the impact of depletion as well as of a 
variation of the external FUV field. 

Results. We find it difficult to simultaneously model the observed '^CO and ''CO emission. The '-CO and [Ci] emission can be explained 
by a PDR model with a external FUV field of 1-0.75 ;fo, but this model fails to reproduce the observed ''CO by a factor of ~2. Adding PAHs 
to the chemical network increases the [C i] emission by 50% in our model but makes [C n] very faint. CO depletion only slightly reduces the 
'-CO and "CO line intensity (by <10% and <20%, respectively). Predictions for the [Ci] ^Pi/i-^^Pi/i, [Ci] 'P^-i^'Pi and '-CO J = 5 ^ 4 
and 4— »3 transitions are presented. This allows a test of our model with future observations (APEX, NANTEN2, HERSCHEL, SOFIA). 
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1. Introduction 

The understanding of Photon Dominated Regions (PDRs) is 
of great interest as they account for much of the millimitre, 
sub-millimitre and far-infrared line and continuum radiation 
from the molecular interstellar medium (ISM) of nearby star 
forming regions to clouds in the diffuse Galactic radiation 
field. Many numerical codes have been developed in order to 
model their emission and to und erstand the physical and chem - 
ical processes governing them (iHollenbach & Tielens Il999h . 
Typically, PDR models use a physical structure (morphology, 
density profile) to calculate a self-consistent solution for the 
chemistry (abundance of species) and the energy balance (heat- 
ing and cooling). Practically, the numerical codes differ in the 
scope of the chemical networks and their reaction rates as well 
as the degree of sophistication with which the heating, cool- 
ing and the radiative transfer are calculated. In many cases 
the goal of the PDR modelling is to constrain the physical 
properties of the line-emitting region (for example, its den- 
sity and column density) and the strength of the external far- 
ultraviolet (FUV) field. PDR models successfully account for 
much of the millimitre, sub-millimitre and far-infrared line and 



continuum emission from Galacti c and extra-galactic PDRs 
(IHollenbach & Tielens|[l997[ll999h . 



The coupling of the energy balance (heating & cooling) 
and the chemistry makes PDR modelling not a straightforward 
task, however In particular, a large set of non-linear equa- 
tions needs to be solved in order to obtain abundance pro- 
files of species. Small changes in the initial elemental abun- 
dances or the reaction rates of individual reactions within their 
sometimes large error margins can result in lar ge variations 



of the abundance profiles of individual species dRolhg et al 
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2007h . In some regions of the parameter space, more than 
one steady-stat e solution to the chemi cal network may exist 
(bistability, e.g. iLe Bourlot et al.lll993b . Additional complica- 
tions arise from intrinsic variations of published PDR mod- 
els which are noted even when identical source parameters 
and chemical react ion rates are used. The benchmarking study 
(IRoIIig et al.ll2007h has shown that even for constant density, 
plane-parallel models, the emerging line intensity of different 
numerical codes vary by up to an order of magnitude for some 
of the fine structure transitions. These discrepancies likely re- 
flect the different numerical methods which are used to solve 
the chemical structure as well as the numerical implementation 
of physical processes. The calibration of PDR models cannot 
be easily done since the observed clouds generally have more 
complex morphologies than assumed in most PDR models. It 
is therefore desirable to study clouds with a well-known geom- 
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Table 2. Standard Model Parameters 



Radius (arcsec) 



Fig. 1. Extinction profile for Barnard 68. Filled dots repre- 
sent the measured Ay calculated by Alv es et al. (2001), ex- 
cluding the south-east prominence seen in the visual extinc- 
tion map. The solid line gives the line-of-sight extinction pro- 
file versus radius for the power-law density profile described in 
Section l4n The dashed line is the same column density profile 
for a model with the same density power-law index but a 10% 
larger radius. 



etry (density profile) which can be used to calibrate PDR mod- 
els. For example, numerous PDR models have been made for 
the Orion Bar because of its approxima t e plane-parallel geom- 
etry ( e.g. Tielens & Hollenbachl 1985 : Hollenbach & Tielens 
1999h . 

Barnard 68 is an example w here the density pr ofile is con- 
strained to a very high precision dAlves et alfcOOlb . and its ge- 
ometry is well matched to a spherically-symmetric PDR model. 
With no young stars nearby, the external FUV field is low, 
of the order of the diffuse Galactic radiation field. Due to its 
proximity to the Op hiuchus complex, w e adopt a distance to 



the cloud of 125 pc dde Geus et al, 



extinction techniques, lAlves et al 



19891) Using near-infrared 



(1200 Ih derived an extinc- 



tion map that was used to constrain the radial density profile 
of the cloud, suggesting that Barnard 68 is consistent with 
a pressure-confined self-gravitating cloud near equilibrium, a 
Bonnor-Ebert sphere. Avoiding the uncertainly concerning the 
geometry and density structure we are able to study in detail 
the impact of the chemistry in the predicted line emission of 
the model. 

Here, we present observations of the three lowest rotational 
transitions of '^CO, the [C i] ^Pi — >^P() fine structure transition 
of neutral carbon, ^^CO 7 = 2 ^ 1, and '^CO 7 = 3 2 and 
compare them to t he KOSMA-t spher ical PDR model origi- 



nally published by iStorzer et al.l (11996 ). The observations are 



summarized in Section |2] Section |3] describes the employed 
model, and the results are discussed in Section]?] In Section ]5] 
the gas heating and cooling is studied and predictions of the in- 
tensity for other important cooling lines are made. A summary 
is given in Section]6] 



Cloud surface density (cm-'') 
Cloud radius (cm) 
Exponent of the density power-law 
Line Doppler width (km i"') 



2.01 X 10^ 
1.875 X 10'^ 
1.96 
0.11 



Table 3. Fractional Abundances" 



He 


0.1 


C 


1.32 X 10-'* 





2.91 X 10"" 


N 


8.5 X 10-5 


"C 


2.08 X 10^*^ 


S 


1.87 X lO^*^ 


Mg 


5.12 X 10-*^ 


Fe 


6.19 X lO^*^ 


Si 


8.21 X 10^' 


PAH'' 


1 X 10-' 



" Abundance relative to H nuclei. 

^ Not included in the reference model. 



2. Observations 

Observations were made for the three lowest rotational transi- 
tions of CO and the [C i] ^Pi — »^Po fine structure transition of 
neutral carbon in Barnard 68. Additionally, we use the '^CO 
J - 3 2 and 7 = 2 — > 1 rotational transitions observe d 
with the CSO telescope and published by B ergin et al. I (l2006h . 
Using these transitions we deliberately limit ourselves to sur- 
face tracers, since the goal of this paper is to model of the PDR 
emission of a cloud in the diffuse Galactic radiation field. 

The high signal-to-noise [Ci] ''Pi^^Po observation was 
made in 2001 with the Submillimiter Wave Astronomy Satellite 
(SWASQ) and has an angular resolution of 4'3. It is pointed 
at a position (Aa - -0^094, A6 =-i-0;2) offset from the col- 
umn density peak of Barnard 68 at a = 17''22 '"38'!6 and 6 - 
-23°49'46'.'0 (J2000) (]Di Francesco et al.]]2002]) and covers es- 
sentially the whole cloud. The r.m.s noise AT^^ is 23 mK per 
velocity channel (Av= 0.623 km s"^). A 3rd-order polyno- 
mial is fitted and subtracted from the spectra. S WAS simultane- 
ously observes the transitions of four species, [C i] ^^Pi — >''Po at 
492 GHz, i^CO 7 = 5 4 at 550.9 GHz, H2O at 556.0 GHz, 
and O9 at 487.2 GHz. No emission is detected for the latter 



three species, however, and .Bergin & Snelll (120021) give the re- 
sulting upper limits on the H2O and O2 column densities. The 
rm.s noise of the '^CO 7 = 5 — > 4 spectrum is 13 mK per 
velocity channel. 

The '^CO 7 = 1 ^ (115 GHz) rotational transition was 
observed in June 2005 with the Mopra 22-m telescope located 
in Australia. A 6' x 6' area was mapped using the on-the-fly 
mode. The data were smoothed to an angular resolution of 1 ' on 
a grid with a 30" spacing. The resulting spectra typically have 
a rm.s noise AT^^ better than 0.5 K per velocity channel (Av= 
0.08 km s"'). A straight line was fitted and subtracted from 
each spectrum. 



See lMelnick et alj ( l2000h for additional details about SWAS. 
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Table 1. Spectral Line Observations Toward Barnard 68 
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Telescope 


Transition 


Frequency 


# of Positions 


Sampling 


%b 


'?mb 


Av 






GHz 


[arcsec] 


[arcsec] 




[kms-'] 




SWAS 


[Ci] 'P, ->^Po 


492.160 


1 


- 


258 


0.9 


0.623 


Mopra 


'-coy = 1 -»o 


115.271 


144 


15 


33 


0.42 


0.08 


KOSMA 


'-CO/ = 2^ 1 


230.538 


49 


60 


130 


0.76 


0.1 


KOSMA 


'-CO J = 2,^2 


345.795 


49 


60 


80 


0.78 


0.294 


KOSMA 


'^CO J = 2,^2 


345.795 


1 




80 


0.78 


0.024 


CSO' 


"CO J = 2^1 


220.398 


65 


24 


33 


0.70 


0.06 


CSO' 


"CO y = 3 ^ 2 


330.587 


144 


24 


22 


0.75 


0.04 



' Published by Bergin et al. (2006") . 



The '^CO 7 = 2-^ 1 (230 GHz) and '^CO 7 = 3^2 
(345 GHz) transitions were mapped simultaneously with the 
230/345 GHz dual-channel receiver at the KOSMA 3-m tele- 
scope toward 49 positions, spaced by 1'. For the '^CO 7 = 2 — > 
1 transition 90% of the spectra have a rm.s noise Ar,^i, smaller 
than 0.3 K per velocity channel (Av= 0.1 km s"'). The '^CO 
7 = 3^2 spectra are smoothed to the 130" angular reso- 
lution of the 230 GHz beam, and 90% of the resulting spec- 
tra have a r.m.s noise AT^n,), smaller than 0.5 K per velocity 
channel (Av= 0.294 km s"'). Typically, a sinusoidal or 3rd- 
order polynomial was fitted and subtracted from each spectrum. 
Because of a small misalignment of both receiver channels the 
pointing of the 345 GHz beam is offset by AAZ=-7" and 
AEL=H-28" in horizontal coordinates relative to the lower fre- 
quency beam. The observations were completed within a rel- 
atively short period of time during the transit of the source. 
Thus, the offset in the horizontal system translates into an off- 
set of approximately Aa - -1" and A5=-28" in equatorial co- 
ordinates, which was corrected in the final map. We estimate 
the accuracy of this correction to be within ~15". The central 
position was also observed in the '^CO 7 = 3^2 line us- 
ing the high-resolution spectrometer (HRS). The calibration of 
the KOSMA data was regularly checked and corrected using 
a source of known intensity (DR21). The observed spectra are 
shown in Figure |2] and a summary of the observations is given 
in Table □ 

3. PDR Model 

We employ th e spherically symme tric PDR model originally 



Table 4. Summary of PDR models 



developed by IStorzer et al] ( 1 19961) (KOSMA-t model). This 



model uses a spherical cloud with a truncated power-law den- 
sity profile and an isotropic FUV radiation field. We adopt a 
density profile of, n{r)=n^(r / r^y for 0.3rc < r < r^, and con- 
stant density, n(r)=«s(0.3)"" in the cloud center (r < 0.3rc), for 
the present paper Here, is the cloud radius and is the den- 
sity at the cloud surface. In our models we adopt a power-law 
exponent of a - \ .96 for the density profile, a cloud surface 
density of 2x10"* cm"^, and a cloud radius of 1.9x10^^ cm. 
This gives a reasonable fit to the measured column density pro- 
file for r < 100", but does not account for the low-column den- 
sity gas close to the surface at 100" < r < 120" (Figure [T]l. 
We tested the response of the emerging CO emission to small 
variations in the model cloud size. An extreme case is shown 
in Figure [U by the dotted line (rc - 2.09 x 10'^ cm), which 



Model 


FUV field 


PAHs 


Depletion 


# 


Xq 






la 


1.0 


N 


N 


lb 


1.0 


N 


Y 


2b 


1.0 


Y 


Y 


3a 


0.75 


Y 


Y 


3b 


0.12 


Y 


Y 



matches the observed Ay at 100" < r < 120", but overes- 
timates the column density for 70" < r < 100". The latter 
model produces CO emission which is typically larger by only 
a few percent, however. The Doppler line-width in the model 
is b-Q.l \ km s"', derived from the observed '^CO 7 = 2 — » 1 
line profile. The gas temperature measur ed using NH^ obser- 
vations toward Barnard 68 is 10 - 16 K ("Hotze l et a"l]|2002a : 



Bourke et al. 1995). This suggests a pure thermal line width 



of 0.08 - 0.1 km s" for CO. The chemical network of the 
model includes H, He, C, O, N, plus a number of heavier el- 
ements, S, Si, Fe, and Mg. '^C is considered in the chemi- 
cal network, including the isotope selective reactions for '^CO. 
Self-shielding of '^CO and '^CO and shielding by H2 against 
ph oto-destruction is implemente d using the shielding factors 
bv Ivan Dishoeck & Black (Il988l) . A summmary of the model 
parameters is shown in Table |2] and the initial fractional abun- 
dances are listed in Table[3] Line intensities and line profiles are 
determined for the [C 11] and the [C i] fine structure transitions 
and the '^CO and '^CO rotational transitions using the radia- 
tive transfer code by Gierens et al. (1992). For a comparison 
to the observed data we smoothed the model intensity distribu- 
tion and line profiles to the angular and velocity resolution of 
our observations. A summary of the models presented in this 
publication is shown in Table H] 



4. Results 

4.1. Reference Model 

In the following, we define a reference model (Model #1 a) with 
a strength of the external FUV field of = ^-^Xo^ w here xq 
is the intensity of the mean interstellar radiation field dDraine 
1978b . Our modelling goal is to match the observed line- 
integrated intensity and line profiles towards the cloud center as 
well as the azimuthally averaged intensity distribution. Table |5] 



4 
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1^1^ Mtel 





Aa n 



-5 5 10 

Velocity (Km/sec) 

Fig. 2. U/9/?er /e/f; Map of the '^CO 7=1^0 emission of Barnard 68. The map is shown in its original resolution of 33". 
The contours run from 40% to 90% of the peak intensity (72.6 K km s"') in steps of 10%. Upper right: Spectra of the ^^CO 
7 = 2-^1 (thick lines) and 7 = 3^2 (thin lines) transitions with a resolution of 130". The velocity range is v = ... 8 km s \ 
and the main-beam temperature shown is Tmb - - 1 ... 5 K . The offset is given in arcminutes and are relative to the peak of the 
dust emission at a = 17''22"'38^6 and 8 = -23°49'46'.'0 (J2000). Lower left: [C i] ^^Pi ^^Pq emission towai'd the central position 
of Barnard 68. The 43 beam size of SWAS covers a large fraction of the cloud. 



shows the [Ci], '^CO, and '^CO integrated intensity towards 
the cloud center (ratio of the model simulation to the observa- 
tions). 

Model #la reproduces the line-integrated intensity toward 
the cloud center for the [Ci] ^Pi ^^Pq, '^CO 7=1^0, 
and 7 = 2 — > 1 transition within 20%. However, the '^CO 
7 = 2 — > 1 and 7 = 3-^2 model line profiles show self- 
absorption which is not observed (Figure[3]l. 

We note that the '^CO 7 = 3^2 mapping observations for 
r < 100" (Figure |6]l have a line intensity larger than the obser- 
vation in Figure [3] made toward the center position. The '^CO 
observations made with a high velocity resolution of ~ 0. 1 km 
s"' show a red line wing, with emission at velocities between 
3.7 and 4.3 km s"' (e.g. left panel of Figure O. This compo- 
nent cannot be modeled in the framework of a single, spheri- 
cal clump; for a comparison to the PDR model simulations we 



therefore calculated the observed line-integrated intensity ex- 
cluding this red line wing. This was not possible for the '^CO 
7 = 3 — > 2 mapping data because of the insufficient velocity 
resolution of the variable resolution spectrometer (VRS), and 
we expect that the line-integrated intensity of '^CO 7 = 3^2 
in Figure |6] is systematically overestimated. We estimate the 
magnitude of this effect to be of the order of 20%, based on 
our velocity-resolved '^CO observations. Even with the line 
wing subtracted, however, the integrated intensity of the '^CO 
7 = 3 — » 2 HRS spectrum is about 40% smaller than that for 
the corresponding VRS observations made at a different epoch. 
The integrated intensities of both observations only marginally 
agree within 3cr, given the error bars of ~16%. Despite of a 
careful check of the pointing and calibration of the observations 
we have to attribute a possible systematic error of up to ~40% 
to the 7 = 3^2 data. This is likely to be due to the low eleva- 
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Fig. 3. Observed '^CO and '^CO 7 = 3 — » 2 line profiles. The reference model (Model #la) and a model which considers the 
effects of depletion (Model #lb) are also shown. 



tion of the source during the observations (Barnard 68 transits 
at < 20°elevation at the location of the KOSMA observatory). 
While the '^CO J - 3 ^ 2 VRS observations are consistent 
with the model, the HRS spectrum appears to be too weak by a 
factor of 0.5. 

The largest discrepancy between the model and the obser- 
vations is noted for the '^CO emission, however, with line- 
integrated intensities in the model being larger by 1.9 and 3.2 
for the y = 2 — » 1 and 7 = 3 — > 2 transitions, respectively. 
Since the density profile of Barnard 68 is well constrained, the 
discrepancies with the observations must have their origin in 
the chemical network or the strength of the FUV radiation field. 

4.2. Depletion 



The observed line emission is significantly overestimated by 
the model for those transitions which probe somewhat deeper 
layers into the cloud surface, most notably the '^CO transi- 
tions. Observations of C'^O and C'^O in Barnard 68 reveal that 
CO, among other molecules, is depleted f rom the gas phase at 
Ay > 5 due to freezing on dust grains ([Hgtzel et al 1 l20Q2bt 
Bergin et al. 20061) . In the present 



Di Francesco et al. 200: 



model, we study the impact of the CO depletion on the line 
emission by truncating the '^CO and '^CO abundance profiles 
at Ay = 5 (setting the CO abundance to zero for Ay > 5; Model 
#lb). This is only a rough approximation, however, it allows us 
to assess the impact of depletion on the line emission. Figure[3] 
shows the '^CO and '^CO 7 = 3 -> 2 line profiles for Model 
#la and Model #lb; the resulting line-integrated intensities for 
the observed transitions are summarized in Table |5] The line- 
integrated intensities are smaller by 10-20% in the model with 
depletion, mainly because the line profiles are narrower. We 
confirm that depletion plays a neglig ible role for low-7 '^CQ 



emission and note a larger impact on both ^^CO transitions 
(20%) and the '^CO 7 = 3^2 transition (5%). However, 
the effect of depletion is much smaller than the factor of 2 to 
3 discrepancy noted between the observations and the model 
line intensities for '^CO. Thus, depletion of CO cannot solely 
account for the low integrated intensity observed for the latter 
transitions. In the following models (Model #2b, #3a, and #3b) 
we account for depletion in the same way as done in Model 
#lb. 

4.3. Polycyclic Aromatic l-lydrocarbons 

Polycyclic Aromatic Hydrocarbons (PAHs) are known to play 
an important role in the carbon chemistry, increasing the 
C column de nsity and decreasing the CO abundance in th e 
cloud surf ace dLepp & Dalgarno 1988l: lBakes & Tielensll998l) . 
Following Ka ufman et al., (J 999.) we included PAHs with an 
initial abunda nce of I x I0~ ^ relative to H nuclei and assuming 
a MRN dMathis etal.lll977b size distribution for PAHs contain- 
ing between 30 and 1500 carbon atoms. We consider neutral 
PAHs in our chemical network, their singly charged variants 
(PAH^, PAH ), as well as PAHs with carbon and hydrogen ions 
adsorbed on their surface (PAHC^, PAHH^). We use a list of re- 
actions provided by M. Kaufman (private com munication) and 
calculate t he reaction rates from th e work of iDraine & Sutin 
(1987,) and lBakes & TielensI (Il994 . A comparison of the C+, 
C°, and CO abundance profiles and the gas temperature in the 
model with and without PAHs is shown in Figure |4] We note 
that the results are relatively insensitive to the initial PAH abun- 
dance as long as PAHs are included. Models with a factor of 4 
larger PAH abundance give similar line-integrated intensities 
(to within 10 %). One of the main impacts of PAHs is that the 
C^ layer at the cloud surface becomes thinner, increasing the 
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Model #1b 
Model #2b 
Model #3a 
Model #3b 



Fig. 4. Left panel The abundance profile of C^, C°, and CO for Model #lb (thin lines) and Model #2b (thick lines) are compared. 
Central panel The abundance profile of CO and '^CO for Model #lb (thin lines) and Model #2b (thick lines). Right panel Gas 
temperature at the cloud surface as a function of Av for models #lb, #2b,#3a, and #3b. 



Table 5. Model comparison with observations toward the Barnard 68 center 



Species Transition 


Telescope 


Resolution (") 


^obs 


Model #13* 


Model #lb 


-^model / ^obs 

Model #2b 


Model #3a 


Model #3b 


[Ci] 'Pi^-^Po 


SWAS 


258 


0.93±5% 


0.97 


0.98 


1.55 


1.39 


0.41 


'^co y = 1 ^ 


Mopra 


60 


0.11 ±9% 


1.08 


0.99 


0.87 


0.82 


0.51 


i^CO 7 = 2^1 


KOSMA 


130 


0.47±2% 


1.24 


1.18 


1.10 


1.00 


0.54 


'^CO 7 = 3^2 


KOSMA/HRS* 


80 


0.72± 15% 


2.26 


2.14 


2.03 


1.56 


0.68 


'-CO i = 3 ^ 2 


KOSMAAVRS* 


80 


1.41 ± 16% 


1.15 


1.09 


1.03 


0.80 


0.35 


'^coy = 2 -> 1 


CSO 


33 


0.19±5% 


1.92 


1.56 


1.38 


1.30 


0.94 


'^^CO 7 = 3-^2 


CSO 


22 


0.25±8% 


3.19 


2.57 


2.13 


1.95 


1.27 


" The intensities are in units of 10 ' ere 


-1 -1 -1 
s cm sr . 















' Observations made at different epochs with the high-resolution spectrometer (HRS) and the variable-resolution spectrometer (VRS). 



total atomic carbon column density. Close to the surface of the 
cloud PAH" neutralizes with C^, while CO is formed some- 
what deeper into the cloud (Ay > 0.4). The larger neutral car- 
bon column density (at Ay < 0.5; Figure |4| and the higher gas 
temperature (Ay < 0.2) gives a stronger [C i] emission com- 
pared to the model without PAHs. The CO line intensity is 
reduced in Model #2b because the emission traces a slightly 
colder region somewhat deeper into the cloud. The lower CO 
abundance at Ay < 0.4 in Model #2b also removes the self- 
absorption in the '^CO 7 = 2 — » 1 line toward the cloud center 
and reduces it for the '^CO 7 = 3^2 line. 

Model #2b gives [C i]/CO line ratios consistent with the ob- 
servations and removes the '^CO 7 = 2 — > 1 self-absorption 
noted in Model #lb. The '^CO line intensities are significantly 
reduced in the model with PAH but are still a factor of 1 .5 - 2.5 
larger than the observations. The [Ci] intensity is also overes- 
timated by a factor of 1.5. 

4.4. Variations of the Radiation Field 

Generally, the line emission in the PDR model is reduced if a 
weaker external FUV radiation field is assumed. We explore the 
impact of lo wer external radiatio n field intensity as a next step. 
According to Bergin et al.l (120061) . the low '-^CO line intensities 
suggest that Barnard 68 is exposed to a FUV radiation field of 



rather low intensity, corresponding to x - O.llxo (0.2 in units 
of the Habing (1968) field). We run two additional models, the 
first with a somewhat smaller x - 0.75a^o (Model #3a) and a 
second with an intensity following Bergin et akl (120061) (x — 
0.12;^'o; Model #3b). All other model parameters are the same 
as for Model #2b. The observed CO line profiles and the model 
line profiles are shown in Figure|5] and the radial profiles of the 
line-integrated intensity are given in Figure |6] A comparison 
between the observed intensity profiles and the model is shown 
in Table [T] In order to quantify their difference, we calculate 
the average ratio of the model simulation to the observations 
for each observed position within r < 100". 

While the line intensities of both species, [C i] and CO, are 
reduced for a smaller the [C i] emission is more sensitive 
to the intensity of the external FUV field than the CO rota- 
tional transitions. A weaker FUV radiation field significantly 
decreases the gas temperature at the cloud surface (Figure |4^) 
but does not change the C° density appreciably (Figure |4j)). 
Generally, the gas temperature in the models is below 23.6 K, 
the excitation energy of the carbon ^Pi level above the ground 
state. In this domain the [C i] line intensity is very sensitive to 
the temperature and thus to x- In contrast, reducing the FUV in- 
tensity results only in small changes of the gas temperature in 
the region where CO is the most abundant carbon species and 
thus has a correspondingly small impact on the optically thick 
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Fig. 5. Observed line profiles compared with Model #lb (dashed), #2b (solid),#3a (thick dotted), and #3b (thin dotted). The 
transitions are indicated in the upper left corner of each panel. 

Table 6. Predicted line-integrated intensities toward Barnard 68 center 



Species Transition 


Resolution - Instrument 


Intensity (10 ' erj 
Model #Ia Model #lb 


-1 -1 
I s cm sr 

Model #2b 


Model #3a 


[Cll] ^P3/2-»'P,/2 


18"-GREAT-S0FIA 


9.49 


9.49 


4.77 


2.44 




13"-HIFI-Herschel 


9.47 


9.47 


4.76 


2.43 


[Cl] 'Pj Pi 


25"-NANTEN2 


2.83 


2.86 


6.65 


5.43 




25"-HIFI-Herschel 












8" -APEX 


2.83 


2.86 


6.62 


5.40 


i^co y = 5 ^ 4 


39"-HIFI-Herschel 


2.49 


2.38 


1.60 


1.37 




56"-CASIMIR-SOFIA 


2.51 


2.39 


1.61 


1.37 


i^co y = 4 ^ 3 


13"-APEX 


2.22 


2.12 


1.83 


1.56 




45"-NANTEN2 


2.23 


2.13 


1.84 


1.56 



CO lines. The exception is Model # 3b with an extremely low 
FUV radiation field of;^' = 0.12;^fo- In this case the cosmic -ray 
heating prevails already at Ay > 0.4, with even lower tempera- 
tures closer to the surface (Figure|4|i. 

The strength of the FUV radiation field also influences the 
self-absorption in the model '^CO J - 2 1 and 7 = 3^2 
line profiles because the latter depends on the gas temperature 
at the cloud surface. As these transitions are sub-thermally ex- 
cited at the lower-density gas near the cloud edge, the reduction 
of the gas temperature in the low-FUV model results in stronger 
self-absorption (Figure|5l thin dotted lines). 



The lower FUV intensity significantly reduces the [Ci] 
and '^CO emission. The model with;^' = O.llxo (Model #3b) 
reproduces the observed '^CO intensities and line-integrated 
intensity profiles to within ~40%, confirming the results by 
Bergin et al. However, this model cannot account for 

the observed [Ci] and '^CO emission, which is underestimated 
by a factor of >2. Additionally, we note that models with a low 
FUV radiation field produce clear self-absorption in the '^CO 
7 = 2^1 and 7 = 3^2 line profiles which is not observed. 
The model with a slightly reduced FUV radiation field (x = 
0.75;^fo; Model #3a) reduces the [Ci] and '^CO emission by 
10% compared to thex - ^-Oxo model. While the '^CO model 
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Fig. 6. Line-integrated intensities as a function of cloud radius. The transitions are indicated in the upper right corner of each 
panel. The observations (error bars) are compared with predictions for Model #lb (dashed), #2b (solid),#3a (thick dotted), and 
#3b (thin dotted). Horizontal error bars in the '^CO 7 = 2 — > 1 and 7 = 3^2 emission correspond to the uncertainly in 
the correction for the receiver beam offset for the KOSMA observations (see Section |2]i. We excluded the lower half of the 
Barnard 68 data from comparison to the CO data because of the significant deviation from the spherical symmetry, visible in the 
dust extinction maps (lAlves et al ] |200lh . 

Table 7. Model comparison with observed cloud intensity profiles 



Species Transition 


Model #lb 


Model #2b 


-^niodel/^ojji; 

Model #3a 


Model #3b 


'^CO 7=1^0 


0.95 


0.82 


0.78 


0.50 


'^CO 7 = 2-^1 


1.10 


0.98 


0.90 


0.48 


'^CO 7 = 3^2 


1.10 


1.00 


0.87 


0.40 


"CO 7 = 2^1 


1.78 


1.49 


1.40 


0.98 


"CO 7 = 3^2 


3.28 


2.60 


2.39 


1.50 



" Sum over all positions with r < 100". 



intensity still marginally agrees with the data (within < 30%), 
the '^CO emission is significantly larger than observed. Any 
further reduction brings the '^CO model intensity closer to the 
observations, but at the same time gives [Ci] and '^CO intensi- 
ties that are too small. 

We notice that the model predictions for the '^CO 7 = 1 — » 
integrated intensity profile show larger discrepancies at radii 
between 90" and 150". This is possibly due to deviations from 
spherical symmetry at angular scales < 60" (note, for example, 
the deviations of the iso-intensity contours from circular sym- 
metry in Figure 121 recall that we excluded the southern half of 
Barnard 68 for the same reason, see Fig. |6]l. These deviations 
are therefore more noticeable in the high-angular resolution 



'^CO 7 = 1 — > map than in the KOSMA maps (angular res- 
olution 130"). We checked for any possible contri bution from 



stray- emission from the Mopra error beam (~60", iLadd et al 



2005h and concluded that this cannot account for a significant 



fraction of the excess emission observed at 90" < r < 150". 
Also note that emission from a low- Ay region at 100 < r < 120 
(Figure [T]) cannot account for this excess (see discussion in 
SectionlO)- 
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Fig. 7. Dominant heating and cooling processes for Model #3a (x - 0.75xo)- The heating and cooling rates are in units of erg s 

cm-'' . 



Table 8. Line-integrated intensity averaged over the projected 
clump area for Model #3a (x = Q.75xo)- 



Species Transition 


Intensity 




(IQ-^ erg S-' cm-- sr-') 


[Cll] 'P3/2^'Pl/2 


3.94 


[CI]3P, ->3po 


3.76 


[Cil^Pj Pi 


6.95 


'^CO J = 3 ^2 


1.19 


'-CO y = 4 -> 3 


1.38 


'-co y = 5 -> 4 


1.13 


'^co y = 2 ^ 1 


0.23 


'^co y = 3 ^ 2 


0.45 


'^co y = 4 3 


0.30 



5. Discussion 

5.1. Influence of the radiation field 

We find that for the given density profile of Barnard 68 it is not 
possible to find a PDR model that matches all observations. The 
main difficulty is to simultaneously match the '^CO and '^CO 
transitions. The weak '■'CO lines suggest a low FUV intensity 
of;^f = 0.12^^0, but underestimate the [Ci] and '^CO intensity 
by up to a factor of ~2 and produce deep self-absorbed '^CO 
line profiles towards the cloud center. The '^CO and [C i] inten- 
sities suggest a moderately reduced FUV field = 1 - 0.75;^o- 
However, this overestimates the '^CO intensity by a factor of 
1.5 - 2.5. A strength of the FUV field of the order of the 
mean interstellar radiation field have been i ndependently in- 
ferred from observations at Ijum and 90/im by Galli et al]|2002 
(~2.5 Go, corresponding to ~ l.Sxo), arguing against a substan- 
tial reduction of the FUV field at the location of Barnard 68. A 
smaller impact is noticed by including PAHs in the chemical 



network of the model. The [C i] intensity is increased by 60%, 
while the CO intensity is reduced by 20%. PAHs help to re- 
move or reduce the self-absorption and make the line profiles 
more consistent with the observations. Finally, we note that de- 
pletion affects the observed transitions only to a small degree. 

Observations of higher-/ CO, [Ci] ^Pa^^Pi, and [Cii] 
transitions are useful to discriminate between the impact of rel- 
evant model parameters and whether the low '^CO intensities 
indeed result from a low FUV field. They will also allow to bet- 
ter assess the impact of PAHs and/or depletion due to freezing 
on dust even at the cloud surface. Predictions for Model #lb, 
#2b, #3a, and #3b are listed in Table |6] for observations with 
present and future observatories. The impact of PAHs is a re- 
duction of the [C ii] intensity by a factor of 2, while the [C i] 
•'P2^"'Pi intensity is increased by the same factor. Because of 
their high excitation energy, both the [Cii] and [Ci] ^P2— >^Pi 
transitions are very sensitive to the FUV field. They are re- 
duced by up to a factor of 0.5 when the FUV field is reduced 
to X = 0.75;i'o, and by as much as a factor of 1.4x10-* for 
X = 0.12;t'o- The '^CO 7 = 5^4 and 7 = 4^3 line inten- 
sities are more moderately reduced. Together, the inclusion of 
PAHs and a somewhat lower FUV field (x = 0.75;t'o) reduce 
them by 40%. 

The detection of the '^CO 7 = 4^3 (APEX, NANTEN2) 
and '^CO 7 = 5-^4 (HlFl Herschel, GREAT SOFIA) transi- 
tions can be obtained within a few minutes or less of observ- 
ing time. A substantially longer integration time is required 
for the [Cii] ^P3/2— >^Pi/2 transitions. The detection of the [Ci] 
^P2— >"'Pi transition might also prove very difficult from ground 
based observatories. However, with HlFl-Herschel, the [Ci] 
^P2— ^■'Pi transition can be detected in less than a minute. 
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5.2. Heating and cooling 

In Figure |7] we show the heating and cooling for different pro- 
cesses as a function of the visual extinction Ay (Model #3a, 
X = O.lSxo)- Photoelectric heating dominates the gas heating 
at Av < 0.8, being replaced by cosmic-ray heating at larger 
depths. The gas cooling is governed by the [Ci] ^Pi^-^Po and 
^P2^^Pi fine structure transitions for Ay < 0.2, and by the 
'^CO (and '^CO) rotational transitions for 0.2 < A^ < 1.2. Gas- 
grain coupling dominates at Ay > 1.2. It is interesting to note 
that the [C ii] transition does not dominate the cooling in this 
case, even very close to the cloud surface (C^ dominates only 
at Av <0.05). 

Table [8] lists the line emission averaged over the projected 
cross section of the clump for Model #3a (x = 0.75;t^o). The 
[Ci] fine structure transitions dominate with a total intensity 
of 1.1x10"^ erg s ' cm"^ sr followed by the '-CO and 
'^CO rotational transitions with a total intensity of 4.7x10"^ 



erg s ' cm ^ sr The [Cii] "^P 



3/2" 



1/2 



fine structure tran- 



sition contributes with a total intensity of 3.9x 10"^ erg s"' 
cm"^ sr"'. This relatively low value is attributed to the low 
FUV radiation field of Model #3a and to the presence of PAHs 
in the chemical network of the model, which significantly de- 
creases the abundance of close to the surface of the cloud. 



6. Summary and Conclusions 

We present simulations with a spherical PDR code to model 
the line emission in Barnard 68. Its approximately spherical 
geometry makes Barnard 68 an ideal test for such a spherically- 
symmetric PDR model. We compare the model results with the 
line profiles and the azimuthally averaged intensity profiles of 
the three lowest rotational transitions of '^CO, '^CO 7 = 3 — » 
2, 7 = 2 — > 1, and the spatially and spectrally unresolved [Ci] 
^Pi— >^Po transition. This is the first time that a spherical PDR 
code was used to model the line emission in a low-FUV cloud 
in such a detail. 

We have tested the impact of PAHs in the chemical net- 
work of the model, the FUV radiation field intensity, and con- 
sidered the eff'ects of depletion. The impact on the emerging 
'^CO, '^CO, [Ci], and [Cii] line emission is quantified. We 
find evidence that PAHs play a role in the chemical network. 
This results in a larger C" layer, while the C^ layer is very thin, 
making the [C ii] emission very faint and diflicult to discrimi- 
nate from emission from the ambient Warm Ionized Medium 
(WIM) in spectrally unresolved observations. Depletion is not 
likely to play an important role in the observed transitions. 

We find it diflicult to simultaneously model the observed 
'^CO and '^'CO emission, with residual differences up to a fac- 
tor of 2 between PDR models and observed line intensities. 
The weak '^CO line intensity seems to require a steep drop in 
either the excitation conditions or in the abundances in a layer 
between the regions which dominate the '^CO and '^CO emis- 
sion, res£ectively.A_model with a low external FUV field (see 
also iBergin et al. I l2006h is clearly incompatible with the '^CO 
and [Ci] observations. In order to reduce the '^CO line emis- 
sion, a shift of the '^CO abundance to a region where cosmic 
rays dominate the gas heating (Av > 0.8) is required. However, 



isotope-selective reactions produce the opposite eff'ect and a 
'^CO abundance corresponding to the '^C isotope abundance 
is reached in a region which is even closer to the cloud surface 
than for '-CO (Fig. IHi. Consequently, the low-7 '^CO emis- 
sion is optically thick at Av = 0.3 - 0.5 where the gas heating 
is still dominated by the external FUV radiation field (Fig.|7]l. 
Including PAHs in the model reduces the gas temperature in 
the region where '^CO is emitted, but this effect turns out to 
be insufficient. The '^CO intensity could be significantly re- 
duced if the gas-grain collisions dominate the cooling closer 
to the surface and the dust temperature was close to ~10K. 
However, a confirmation of the low observed '^CO intensity 
is also highly desirable. Together with observations of mid-7 
'^CO, [Ci], and [Cii] transitions with present and future tele- 
scopes (NANTEN2, APEX, SOFIA, Herschel), this will allow 
a test of our results and assess the importance of different model 
parameters. 
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